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^ I We analyze the topology, lifetime, and emissions of a torus around a black 

fS| ! hole formed in hypernovae and black hole-neutron star coalescence. The 



torus is ab initio uniformly magnetized, represented by two counter oriented 
current-rings, and develops a state of suspended accretion against a "mag- 
^ _ netic wall" around the black hole. Magnetic stability of the torus gives 

I rise to a new fundamental limit S^/Sk < 0.1 for the ratio of poloidal mag- 

^ I netic field energy-to-kinetic energy, corresponding to a maximum magnetic 

field strength 5, ~ lO^^G (TMq/Mh) (eM^/i?)^ (Mt/0.03Mh)^^^ The life- 
I time of rapid spin of the black hole is effectively defined by the timescale of 

^ [ dissipation of black hole-spin energy Erot in the horizon, and satisfies T ^ 

40s{MH/7MQ){R/6MH)^{0mMH/MT) for a black hole of mass Mh surrounded 
by a torus of mass Mt and radius R. Erot of the black hole. The torus con- 
verts a major fraction Eg^/Erot ~ 10% into gravitational radiation through 
a finite number of multipole mass-moments, and a smaller fraction into MeV 
neutrinos and baryon-rich winds. At a source distance of lOOMpc, these emis- 
sions over = 2 X 10^ periods give rise to a characteristic strain amplitude 
VNhchar — 6 X 10~^^. We argue that torus winds create an open magnetic 
flux-tube on the black hole, which carries a minor fraction Ej/Erot — 10^^ in 
baryon-poor outflows to inflnity. We conjecture that these are not high-sigma 
outflows owing, in part, to magnetic reconnection in surrounding current sheets. 
The fraction Ej/Erot ~ {^/^){Mh/ R)"^ is standard for a universal horizon half- 
opening angle 6h — Mh/R of the open flux-tube. We identify this baryon poor 
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output of tens of seconds with GRBs with contemporaneous and strongly cor- 
related emissions in gravitational radiation, conceivably at multiple frequencies. 
Ultimately, this leaves a black hole binary surrounded by a supernova remnant. 



1. Introduction 

Black holes surrounded by a magnetized torus or disk are believed to constitute the 
central engines that power various high-energy sources, notably active galactic nuclei, galactic 
microquasars, and gamma-ray bursts. The latter systems are thought to be the outcome 
of catastrophic events such as core-collapse in massive stars and black hole-neutron star 
coalescence (Eichler, et al. 1989; Woosley 1993; Paczynski 1991; Paczyhski 1998), and are of 
interest as potentially the most extreme and short-lived black hole-torus systems. 

We describe a theory for the topology, hfetime and emissions of a black hole-torus sytems 
as a function of three parameters: the mass Mh of an extreme Kerr black hole, the radius R 
and the mass Mj- of the torus. We shall do so largely by studying the torus by equivalence to 
pulsars, in both topology and ms rotation periods. The energy emissions arc powered by the 
spin-energy of a Kerr black hole (Kerr 1963). Most of the black hole- luminosity - the rate at 
which the black hole deposits energy into its surroundings in all channels - is incident on the 
torus, which hereby creates a major energy output of the system. A minor energy output is 
released in baryon-poor outflows through an open magnetic flux-tube along the spin-axis of 
the black hole. The hfe-time of these black hole-torus systems is identified with the hfetime 
of rapid spin of the black hole in a state of suspended accretion (van Putten & Ostriker 2001) 
against a "magnetic wall" around the black hole (van Putten 1999). The suspended accretion 
state results from a strong coupling of the torus to the spin-energy of the black hole. We 
point out that this mechanism is based on a uniform magnetization of the torus, represented 
by two oppositely oriented current rings (van Putten 1999). This magnetization is a natural 
outcome of both black hole-neutron star coalescence and core-collapse in hypernovae. 

In this paper, we quantify (1) the lifetime of rapid spin of the black hole in terms of 
a new magnetic stability criterion for the torus; (2) baryon-rich outflows from the torus at 
MeV temperatures; and (3) the fraction of black hole spin-energy in baryon-poor outflows 
through an open magnetic flux-tube on the black hole, created from outer layers of the inner 
torus magnetosphere by these torus winds. 

The torus is luminous in various channels, which are strongly correlated by the prop- 
erties of the torus: in gravitational radiation, winds, thermal and MeV neutrino emissions 
(van Putten 2001b; van Putten & Levinson 2002). Calorimetric constraints on the torus 
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winds hereby obtains predictions for the proposed emissions in gravitational radiation, while 
calorimetry on the gravitational wave emissions by upcoming gravitational wave-experiments 
obtains a method for identifying Kerr black holes as objects in nature (van Putten & Levinson 
2002). 

Gravitational radiation forms a major output of the system and the dominant output 
of the torus (van Putten 2001b; van Putten & Levinson 2002). This is emitted by a finite 
number of multipole mass moments in a torus of finite slenderness, due to the Papaloizou- 
Pringlc instability (Papaloizou & Pringle 1984; van Putten 2002a) and, conceivably, other 
wave-modes in the torus. These emissions arc candidate sources for the upcoming grav- 
itational wave-experiments by laser interferometric instruments LIGO (Abramovici et al. 
1992), VIRGO (Bradaschia et al. 1992), TAMA (Masaki et al. 2001) and GEO (e.g., Schutz 
& Papa (1999)), or by any of the bar or sphere detectors presently under construction. The 
frequency in gravitational radiation is determined by the keplerian frequency of the torus. 
The latter is strongly correlated to the output energy in torus winds. This suggests the 
possibility of performing calorimetry on the impact of these torus winds on the remnant 
stellar envelope (van Putten 2002b) and on hypernova remnants, in order to constrain the 
expected frequency in gravitational radiation. 

Baryon-poor outflows form a small fraction of the total output from the black hole 
through an open magnetic fiux-tubc (van Putten & Levinson 2002). We here attribute 
the formation of this open fiux-tube to powerful baryon-rich torus winds, driven from the 
surface by escaping McV neutrinos. This neutrino output provides the dominant cooling 
mechanism of the torus during the suspended accretion state, in addition to the energy 
release in gravitational waves. The fraction of rotational energy thus released is proportional 
to 9jj, which is standard for a universal horizon half-opening angle 9h of the inner tube. 
These outflows are probably not high-sigma, owing to magnetic reconnection in the interface 
between the baryon poor and baryon-rich winds in a surrounding outer flux-tube. The 
half-opening angle is possibly related to curvature in poloidal topology (van Putten 2002b) . 

Tentative observational evidence for a black hole-luminosity incident into surrounding 
matter in case of supermassive black holes is found in MCG-6-30-15 (Wilms et al. 2001), and 
in case of stellar mass black hole-candidates in the galactic source XTE J1650-500 (Miller 
et al. 2002). Independent observational evidence of the presence of magnetic fields remains 
elusive. We point out, however, that in our model heating of the torus is due to viscous 
shear between its inner and outer faces, possibly in the form of magnetohydrodynamical 
turbulence, which is different from electric dissipation as envisioned in Wilms et al. (2001). 

In §2 and §3 we describe the topology and stability of the magnetospherc of a torus 
formed in black hole-neutron star coalescence and hypernova, the suspended accretion state 
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around rapidly rotating black holes and the secular time-scale of its rapid spin. The formation 
of a finite number of multipole mass moments by the Papaloizou-Pringle instability in tori 
of finite slenderness is summarized in §4. Energy emissions in various channels by the torus 
are described in §5. We calculate the neutrino driven mass loss rate from the torus' surface 
in §6. These baryon-rich torus winds have some implications for energy extraction from the 
black hole and the creation of open magnetic flux-tubes from outer layers in the inner torus 
magnetosphere. §7- §9 describe the structure of the inner flux-tube supported by the black 
hole, and dissipation by magnetic reconnection in its interface with the surrounding outer 
flux-tube. We conclude with observational consequences in §10. 



2. Formation and structure of the torus magnetosphere 

A torus formed in core-collapse of a massive star or black hole-neutron star coalescence 
will be magnetized with a remnant of the progenitor magnetic fleld. An aligned poloidal 
magnetic fleld in the progenitor star provides a magnetic moment density in the torus, 
aligned with its axis of rotation. Equivalently, the magnetic fleld in the torus is produced 
by two concentric current loops with opposite orientation. In the case of a torus formed 
from the break-up of a neutron star around a black hole, these two current loops form out 
of a single current loop representing the magnetization of the neutron star, upon stretching 
the latter around the black hole followed by a reconnection (see Fig. 1). Conceivably, the 
magnetic field in the torus is amplified by winding or a dynamo process. The stability of 
this magnetic field configuration is discussed in §3, and sets an upper limit of 5 ~ lO^^G on 
the field-strength. 

Below, wc outline the resulting magnetosphere of the black hole torus system in the 
force-free limit. For illustrative purposes we consider first the poloidal topology of the 
vacuum magnetic field configuration. We proceed by discussing the more realistic situation 
of a force-free magnetosphere, and show that the the inner and outer face of the torus are 
each equivalent to a pulsar with, however, generally different angular velocities. 

In the subsequent sections we shall consider some essential energetic aspects in greater 
detail. In particular, it will be shown that appreciable matter outflows are expected along 
open fleld lines to inflnity, and a reconnection boundary is expected to form near the rotation 
axis. This may give rise to a non-force-free magnetic fleld in those regions. 
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Fig. 1. — A uniformly magnetized torus around a black hole (C) is represented by two 
counter-oriented current rings in the equatorial plane. It forms a common end point of both 
core-collapse (A1,B1,C) and black hole-neutron star coalescence (A2,B2,C). Core-collapse 
(Al-Bl) in a magnetized star results in a uniformly magnetized, equatorial annulus (C); tidal 
break-up (A2-B2) wraps the current ring representing the magnetic moment of a neutron 
star around the black hole which, following a reconnection, leaves the same (C). 
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Fig. 2. — The poloidal topology of magnetic flux-surfaces in vacuum is illustrated in flat 
space-time, produced by two counter-oriented current rings representing a uniformly magne- 
tized torus (center). The dashed line is the separatrix between the flux-surfaces supported 
by the inner and the outer faces of the torus. A Kerr black hole develops an equilibrium 
magnetic moment which preserves essentially uniform and maximal magnetic flux through 
its horizon. 
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2.1. Topology of the vacuum magnetic field 

A uniform magnetization of the torus is approximately described by two counter-oriented 
current rings in the equatorial plane. A third current loop is associated with the black hole, 
representing its equilibrium magnetic moment in its lowest energy state (van Putten 2001b). 
This induced magnetic moment is oriented antiparallel to the magnetic moment of the torus, 
facilitating an essentially uniform and maximal horizon flux at arbitrary spin-rates. In flat 
space-time, the magnetic field produced by a superposition of current rings can be calculated 
analytically (Jackson 1975). The topology of this ab initio flat-spacetime vacuo magnetic 
field is shown in Fig. 2. As seen, at large radii (compared with the radius of the outer current 
ring) it quickly approaches a dipole solution. In the inner region the field lines intersect the 
horizon, giving rise to a strong coupling between the black hole and the inner face of the 
torus. This topology and flux-distribution is preserved in the face of general relativistic 
effects, as a result of the equilibrium magnetic moment of the black hole. 

2.2. Equivalence to pulsars 

By vacuum break-down, the flux-surfaces will evolve with electric charges to a largely 
force-free state similar to pulsars (Goldreich & Julian 1969). As a result, a magnetosphere 
develops which consists of conductive flux-surfaces and magnetic winds. The torus hereby 
supports an inner and an outer torus magnetosphere. These are equivalent in poloidal 
topology to pulsar magnetospheres, wherein the horizon of the black hole is equivalent to 
a compactified infinity with non-zero angular velocity (Fig. 3). This equivalence implies 
similar, causal interactions by magnetic winds acting on the inner and outer face of the 
torus. 

In the force-free limit, the fiux-surfaces in the outer/inner torus magnetosphere assume 
rigid corotation with the outer/inner face of the torus, by no-slip boundary conditions on 
the surface. Field-lines in the vicinity of the torus and the equatorial plane form a 'bag' 
of closed field-lines (both on the inner and the outer face) with no-slip/no-slip boundary 
conditions. The last closed field-line of the outer torus magnetosphere reaches the light- 
cylinder associated with the angular velocity of the outer face of the torus, similar to the 
last closed field-lines in pulsar magnetospheres. The last closed field-line of the inner torus 
magnetosphere reaches the inner light-surface (Znajek 1977) associated with the angular 
velocity of the inner face of the torus (Fig. 2 of van Putten (1999)). Beyond, field-hnes are 
open, and extend to infinity or to the horizon of the black hole with no-slip/slip boundary 
conditions. The former are created by torus winds which cross the outer light cylinder. 
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Spin-up Spin-down 

Fig. 3. — Lower left. The inner face of the torus (angular velocity and the black hole 
(angular velocity ^Ih) is equivalent to a pulsar surrounded by infinity with relative angular 
velocity fin — By equivalence in poloidal topology to pulsar magnetospheres, the inner 

face receives energy and angular momentum from the black hole as a causal process, whenever 
Qh — ^+ > 0. Lower right. The outer face of the torus (angular velocity Q ) is equivalent 
to a pulsar with angular velocity and always looses energy and angular momentum, by 
the same equivalence. 
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2.3. Suspended accretion against a magnetic wall 

Most of the black hole-luminosity is incident on the torus (van Putten 1999). The 
torque exerted on the inner face of the torus by the black hole is obtained by integrating the 
angular momentum flux C = F'^^F^g/An over the section of the horizon which is threaded 
by magnetic field lines that are anchored to the torus: 

T+^Ati f ' ^Cde = {Qh -^+) I ' ^sineiF^efde, (1) 

where eq. (A18) has bee used. Here 9h is the angle of the last field line that connect the 
torus and the horizon (see Fig. 3), and Q,h are the angular velocities of the inner face of 
the torus and the black hole, respectively, and the metric components, p, E are defined in 
the appendix. In terms of the net poloidal magnetic flux associated with the open field lines 

in the torus, 2ttA^ we may write = [VLh — Q^)ffjA'^. Formally fn is defined through eq. 
(1). Likewise, the torque exerted on the outer face by field lines that extend to infinity can 
be expressed as: T_ = Q^f^A^, where is the angular velocity of the outer face of the 
torus, and is the fraction of open field lines that extend to infinity. It is seen that when 
Qh > ^+ angular momentum is transferred from the black hole to the torus, tending to spin 
up the inner face, whereas in the slowly rotating case {fin < ^+) the black hole receives 
angular momentum from the torus (Fig. 4). The outer face always loses angular momentum 
via a wind to infinity. By mechanical work, the magnetic torus winds to infinity and into 
the horizon carry outgoing luminosities 

L± = Ji±r±. (2) 

The inner face thus receives a fraction D,+/D,h of the maximum extractable power flHT+ 
(Thorne et al. 1986); the rest dissipates on the horizon. 

The positive and negative torques exerted on the inner and outer faces of the torus in 
the rapidly rotating case, give rise to a differential rotation of the torus. In steady state, a 
flow of angular momentum from the inner to the outer face is accomplished through shear 
forces. The nature of the latter is presumably electromagnetic (Balbus & Hawley 1992), 
stimulated by magneto-hydrodynamic instabilities (see §4) . As a consequence, a fraction of 
the black hole spin down energy that is intercepted by the inner face is dissipated in the 
torus, heating it to MeV temperatures. 

2.4. Frame-dragging and electric charge distribution 

The electric charge distribution can be obtained from Maxwell's equation: F*^^ = 47rj*. 
Assuming the radial magnetic fleld, Br — —F^g/CE sin 9), to be independent of 9, we obtain 



light cylinder 



Fig. 4. — Schematic illustration of the poloidal topology of the magnetosphere of a torus 
surrounding a rapidly rotating black hole. Balance of the input on the inner face and the 
output through the outer face and other losses results in a suspended accretion state. The 
inner torus magnetosphere hereby represents a magnetic wall around the black hole, whereby 
the torus receives energy and angular momentum from the black hole by equivalence to 
pulsars. The associated horizon Maxwell stresses (Blandford & Znajek 1977) and Maxwell 
stresses on the inner face of the torus (van Putten 1999) are mediated by poloidal currents, 
which close over a vacuum gap in an annulus of vanishing magnetic field. The torus hereby 
receives a major fraction of the black hole-spin energy, catalyzing this into gravitational 
radiation, winds, thermal and MeV neutrino emissions. A baryon-poor inner flux-tube serves 
as an artery for a minor fraction of black hole-spin energy. The dashed line indicates the 
light cyhnder of the outer face of the torus. The lifetime of the system is set by the hfetime 
of rapid spin of the black hole. 
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Fig. 5. — Shematic of the charge-distribution in a torus magnetosphere around Kerr black 
holes, showing two bifurcations from non-rotation (dotted line). Frame- dragging introduces 
a sign-change in flux-tubes which rotate slower than the black hole (a). Black holes rotating 
faster than the torus produce an inner light surface, and hence an annulus of -B = (b). 
Equilibrium magnetic moments of the black hole and the torus correspond to charges Qh 
and Orr with nnnnsitf! sip"n (Or^ is analnernu.s tn nul.sar rhaTres Cnhen et a,l. f1 9751V 



- 12 - 



from eq. (A19) the Goldreich- Julian charge density near the horizon, on field lines that 
emanate from the inner face: 



2.i_ + (3)Br cose 



where P is defined below eq. (Al), and equals —^h on the horizon. We find that for 
Qh > the charge density changes sign on field lines threading the horizon. The same 
holds true also for the inner fiux tube that extend to infinity (see discussion following eq. 
[44]). 

The charge-distribution (3) represents a bifurcation from the magnetosphere around 
a Schwarzschild black hole. Ingoing horizon boundary conditions and outgoing boundary 
conditions at infinity may hereby carry a continuous electric current. If the black hole 
rotates slower than the inner face of the torus, the inner light-surface is absent and the sign 
of the charge-distribution near the torus carries through to the horizon (Fig. 5a). If the 
black hole rotates faster than the inner face of the torus, the inner light-surface becomes 
apparent (a bifurcation from slow rotation) and, by forgoing arguments, introduces a sign 
change in the charge-distribution (Fig. 5b). 



3. The lifetime of rapid spin of the black hole 

A magnetized torus of a few tenths of solar masses around a stellar mass black hole of 
about 7Mq is subject to magnetic self-interaction and a stabilizing tidal interaction in the 
central potential well. We here derive limits on the average strength of a poloidal magnetic 
field which can be supported by the torus. We do not consider the problem of stability of 
a poloidal magnetic field itself, such as the magnetorotational (MRI) instability. A hmit on 
the energy in poloidal magnetic field defines a lower bound on the lifetime of rapid spin of 
the black hole. 

In regards to wave-motion within the equatorial plane, the contribution of poloidal 
magnetic fields is that of magnetic pressure, which is generally stabilizing on the motion of 
the fiuid. In regards to poloidal wave-motion, a poloidal magnetic field generally conspires 
towards instabilities. This can be calculated by partitioning the torus in a finite number of 
fiuid elements with current loops, representing local magnetic moments. The two leading- 
order partitions are shown in configurations C and Bl of Fig. 1, for which we derive critical 
magnetic field-strengths. The first is subject to magnetic tilt instability between the inner 
and the outer face, and the second is subject to a magnetic buckling instability. 
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Fig. 6. — A uniformly magnetized torus is in its highest magnetic energy state. The two 
alternative leading-order partitions of the current distribution (top) have unstable poloidal 
modes, described by a relative tilt between the two current rings towards alignment (A) 
or towards buckling (B), characterized by perturbations out of the equatorial plane with 
poloidal angles 9i ^ 9j. We consider vertical displacements of fluid elements along a cylinder 
of radius R. Their wave-modes are stabilized by the tidal gravitational held of the central 
black hole, provided the magnetic fleld-strength remains below a critical value on the order 
of lO^^G. This gives rise to a minimim lifetime of rapid spin of the black hole on the order 
of tens of seconds, consistent with the redshift-corrected durations of long GRBs. 
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3.1. A magnetic tilt instability 

Following C in Fig. 1, consider the magnetic interaction energy of a pair of concentric 
current rings, given by 

Ui,(e) ^ -liB cose. (4) 

Here, is the magnetic moment of the inner ring, B is the magnetic field produced by the 
outer ring, and 6 denotes the angle between /j, and B. Note that U^{9) has a period 2%, is 
maximal (minimal) when /i and B are antiparallel (parallel; see Fig. 6). Consider tilting a 
fluid element of a ring out of the equatorial plane to a height z approximately along a cylinder 
of radius R. (This is different from tilting a rigid ring, whose clenioiits move on a sphere). A 
tilt hereby changes the distance to central black hole to p = \/rF+^ ~ R{1 + z'^/2R^). In 
the approximation of equal mass in the inner and outer face of the torus, simultaneous tilt 
of one ring upwards and the other ring downwards is associated with the potential energy 

MtMh ( , 1.2/ 



with tan(^/2) = zjR., where we averaged over all segments of a ring. Note that Ug{&) has 
period tt and is minimal when ^ = 0. Stability is accomplished provided that the total 
potential energy \J{&) — 17^(9) + Ug{9) satisfies 

The potential U{9) is shown in Fig. 7, which shows the bifurcation at 6 = 1/12 of the stable 
equilibrium ^ = into an unstable equilibrium with the appearance of two neighboring stable 
equilibria at non-zero angles. The bifurcation point is therefore second order. Nevertheless, 
the torus may become nonhnearly unstable at large angles (6 >> b*). We therefore consider 
below the physical parameters at this bifurcation point. 

For two rings of radii R± with {R+-R_)/{R+ + R_) = 0(1), we have C/^ ~ ^B'^R^ cos9, 
so that (6) gives B^M% = {1/A){Mh/R)^{Mt/Mh), or 

The critical value of the ratio of poloidal magnetic energy {Eb — fsB'^R^/d) to kinetic energy 
{MtMh/2R) in the torus becomes 
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Normalized potential energy in tilt 
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Fig. 7. — The potential energy associated with a poloidal tilt angle between the inner 
and the outer rings is the sum of a magnetic moment-magnetic moment interaction plus a 
tidal interaction with the central potential well of the black hole. It is shown for various 
normalized magnetic field-energies h — SBj^k- The equilibrium ^ = becomes unstable 
when (PU{9)/d9'^ < 0, corresponding to a bifurcation into two stable branches of non-zero 
angles beyond h > 1/12. 
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where fs denotes a factor of order unity, representing the volume of the inner torus magneto- 
sphere as a fraction of AtiR^/S. We emphasize that the hmit (8) is fundamental, independent 
of the mass of the black hole, and the mass and radius of the torus. This result can equiv- 
alently be attributed to stable balance of Lorentz forces against tidal forces, preventing a 
small misalignment to cause a tilt. 

The poloidal magnetic field introduces an anisotropic pressure tensor. At the critical 
magnetic field-strength (7), pressure components in the equatorial plane are predominantly 
magnetic rather than thermal at MeV temperatures. The poloidal pressure components 
(along the magnetic field-lines) arc thermal pressures. In equilibrium, these poloidal pressure 
components are unaffected by poloidal tidal forces if the poloidal fiux-surfaces assume a 
spherical shape inside of the torus. 

The magnetic field may be generated in response to the power received from the black 
hole. If so, then dSs/dt < Lh, where (Thornc, et al. 1986; van Putten 1999) 

Lh ^ ^vB'M]j, (9) 

denotes the black hole- luminosity expressed in terms of a fraction r] = Qt/^h of the angular 
velocity of the torus to that of the black hole. The discrepancy Lh — dSs/dt is carried away 
in the various channels as described in the suspended accretion state, and will be zero in 
equilibrium. This gives rise to a minimum e-folding time 




A critical magnetic field of Be ~ lO^^G is hereby reached on a time-scale of at least a few 
seconds. The suspended accretion state may hereby be intermittent on an intermediate time- 
scale of seconds, associated with magnetic field build-up powered by the rotational energy 
of the black hole. 

Most of the rotational energy is dissipated in the horizon, creating entropy S for a black 
hole temperature Th at a maximal rate ThS ~ (Thorne, et al. 1986). The lifetime 

of rapid rotation of the black hole becomes effectively the timescale of dissipation of black 
hole-spin energy Erot — Mh/^ in the horizon, i.e.: 

T.^>40sMm7^V (11) 

where eq. (7) has been employed. 

The suspended accretion state, conceivably intermittent on a timescale (10), lasts for 
a duration similar to (11) until the ISCO reaches the torus or until Qh — ^t, whichever 
comes first. 
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3.2. A magnetic buckling instability 

We partition the magnetization of the torus into N equidistant fluid elements with 
dipole moments, fii = fi/N = {1/2)BR^/N. We consider the vertical degree of freedom, 
of fluid elements which move to a height z above the equatorial plane. By conservation of 
angular momentum, this motion is restricted to a cylinder of constant radius. Their position 
vectors in and off the equatorial plane will be denoted by 

= (i?cos0j, i?sin0j, 0), = {Rcos4>i, RsiiKpi, Zi), (j)i = 2m/N. (12) 

A fluid element i assumes an energy which consists of magnetic moment-magnetic moment 
interactions and the tidal interaction with the central potential well. The total potential 
energy of the i— th fluid element is given by 

u, = y^^^\ ose,j + u^iOi), (13) 

^j\ 

where B' — B/N* denotes the magnetic field-strength of a magnetic dipole at distance d — 
27: R/N, 6ij denotes the angle between the i— th magnetic moment and the local magnetic field 
of the j— th magnetic moment, and Ugi — —{MtMh / RN){l—z1 /2R^) the tidal interaction of 
the i— th fluid element with the black hole. Here, N* is a factor of order N which satisfies the 
normalization condition EjC/j = —^B (in equilibrium). Upon neglecting azimuthal curvature 
in the interaction of neighboring magnetic moments, we have a magnetic moment-magnetic 
moment interaction 

fliB COS % ~ — 1 - 

\ri - Vj\-^ V- 3\ \ 

where cty = {zi — Zj)/d, 

cos% = -^(l-4)/(l + 4)c^ - (1 - 4) , (15) 

and |rj — r^l ~ \i — j\d (l + afj-/2|i — jp) . 

We shall use the small a small amplitude approximation, whereby z.i/R = tan^j ^ 9i. 
We study the stability of this configuration, to derive an upper limit for the magnetic field- 
strength. An upper limit obtains by taking into account only interactions between neigh- 
boring magnetic moments. (The sharpest limit obtains by taking into account interactions 
between one magnetic moment and all its neighbors.) Thus, we have N* — 2N and consider 
the total potential energy 

U, = ^^i-j\=, {l - \a^)j + U,{e,), (16) 



1 + 



2K-j| 



(14) 
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where aij — N{9i — 9j)/2n. The Euler-Lagrange equations of motion are 



MtR. dUi 

This defines the system of equations for the vector x = (6'i, 6'2, ■ • • , On) given by 

/ 2 -1 ••• -1 \ 



(17) 



MtR MtMh 5aB 
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■1 ••• -1 2 



(18) 



The least stable eigenvector is x = (1,— 1,1,---,— 1) (for N even), for which the critical 
value of the magnetic field is 



5 \Mh V R 



(19) 



This condition is very similar to (7), and gives a commensurable estimate 



£b 

and lifetime of rapid spin of the black hole. 



1 

15 



(20) 



We propose to identify the redshift-corrected durations of long gamma-ray bursts and 
their contemporaneous emissions in gravitational radiation (Fig. 8) with this new secular 
time-scale (11) of tens of seconds, set by the stability- limited energy of the poloidal magnetic 
field supported by the torus. 

A high-order approach can be envisioned, in which the inner and outer face of the 
torus are each partitioned by a ring of magnetized fiuid elements. This is of potential 
interest in studying instabilities in response to shear, in view of the relative angular velocity 
— Q_ > 0. Magnetic couphng between the two faces of the torus through aforementioned 
tilt or buckling modes inevitably leads to transport of energy and angular momentum from 
the inner face to the outer face of the torus by the Rayleigh criterion. A detailed discussion 
of this instability falls outside the scope of the present paper. 



4. Multipole mass moments in a torus in suspended accretion 



As explained in §2, the torus is sandwiched between a magnetic wall around the black 
hole and an outer torus magnetosphere. It is hereby subject to competing torques on the 
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Fig. 8. — Shown is the histogram of redshift corrected durations of 27 long bursts with 
individually determined redshifts from their afterglow emissions. We identify these long 
durations with the lifetime of rapid spin of a Kerr black hole in a state of suspended accretion. 
These durations are effectively defined by the rate of dissipation of black hole-spin energy in 
the horizon, which is subject to a new magnetic stability criterion for the torus. (Reprinted 
from M.H.P.M. van Putten, 2002, ApJ, 575, L71). 
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inner and the outer face, which introduces heating and, possibly, turbulent magnetohydro- 
dynamical flow. This is stimulated by a flnite number of multipole mass moments produced 
by the Papaloizou-Pringle instability in a torus of finite width. The net rate of dissipa- 
tion can be calculated in the suspended accretion state, from balance of energy and angular 
momentum in combination with emissions in gravitational radiation and magnetic winds. 

A torus tends to develop instabilities in response to shear, which can be studied analyt- 
ically in the approximation of incompressible fluid about an unperturbed angular velocity 



Q = a, 



(21) 



where q > 3/2 denotes the rotation index and a = {R^ + -R-)/2 denotes the major radius of 
the torus. In the inviscid limit, irrotational modes in response to initially irrotational per- 
turbations to the underlying flow (vortical if g 7^ 2) shows the Papaloizou-Pringle instability 
(Papaloizou & Pringle 1984) to also operate in wide tori (van Putten 2002a). The neutral 
stability curves of the resulting buckling modes are described by a critical rotation index 
Qc — Qc{^,'>tt) a-s a function of the slenderness parameter 5 — b/2a, where b — (i?_ — -R+)/2 
denotes the minor radius of the torus and a — (i?_ -|- -R+)/2 the mean radius (Fig. 9), where 
m denotes the azimuthal wave-number. 

It will be appreciated that the torus is m = stable for perturbations of its radius (in 
the mean). This is due to the frozen-in condition of magnetic flux-surfaces. In terms of 
angular momentum transport, the horizon magnetic flux oc {M/af hereby deflnes a black 
hole-to-torus coupling which is dominant over the coupling oc — M/ a^/^ of the torus to 
inflnity. 

Quadratic flts to the stability curves are 



5c(5, m) 
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(22) 



Instability sets in above these curves, stability below. For m — 2, the critical value Qc — 2 
obtains for 5 = 0.16, associated with the Rayleigh stability criterion for the azimuthally 
symmetric wave mode m — 0. For large m, we use the numerical result of critical values 
5 = 0.28/m for g = 2. 

The central pressure of the torus is produces by differential rotation, due to a state of 
super-keplerian motion of the inner face and a state of sub-keplerian motion of the outer face. 
This pressure is balanced by both magnetic and thermal pressure at MeV temperatures. 
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Fig. 9. — Diagram showing the neutral stabiUty curves for the non-axisymmetric buckhng 
modes in a torus of incompressible fluid, as an extension of the Papaloizou-Pringle instabil- 
ities to flnite slenderness ratios 6 = b/2a > 0, where b and a denote the minor and major 
radius of the torus, respectively. Curves of critical rotation index Qc are labeled with az- 
imuthal quantum numbers m — 1,2,.., where instability sets in above and stability sets 
in below. The central pressure produced by super- and subkeplerian motions of the inner 
and outer face of the torus is balanced by magnetic and thermal pressure at MeV temper- 
atures. At finite slenderness < b/a << 1, this produces instability to a finite number of 
Papaloizou-Pringle modes (Reprinted from M.H.P.M. van Putten, ApJ, 575, L71). 
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which gives rise to finite slenderness 5 > and hence instability to a finite number of 
Papaloizou-Pringle modes. Thermal pressure alone gives the estimate (adapted from van 
Putten (2002a)) 

Magnetic pressure enhances this estimate to larger values of q. This shows that a torus of 
finite slenderness {6 < 0.1) and a radius around 5M is unstable {q > \/3) to the formation 
of a finite number of multipole mass moments by the Papaloizou-Pringle instabihty. A torus 
with multiple mass moments potentially defines a gravitational wave-spectrum consisting of 
several lines. 



5. Energy emissions by the torus 

The suspended accretion state in the case of a symmetric flux-distribution, given by 
equal fractions of open magnetic fiux on the inner and the outer face, gives rise to remarkably 
simple expressions for the predicted energy output in the relevant channels. To leading order 
the angular velocities of the inner and outer faces are given in terms of the mean angular 
velocity of the torus, D.t = + f^+)/2, and the slenderness ratio 5, as Q± = ^t{'^ ± S). 
Equation 2 implies that for a small slenderness ratio, L_ ~ tjL^, where f] = Qt/^h denotes 
the ratio of the angular velocity of the torus to that of the black hole. In the limit of strong 
magnetohydrodynamical viscosity and small slenderness ratio we then have the asymptotic 
resuhs (van Putten 2002b) 

Egu,/Erot^V, K/Erot^v'^, Ed/Erot^V^, (24) 
where Eg^, E^j, are defined below. 



Gravitational radiation. The major energy output from the torus is in gravitational radiation, 

^,.^6xlO^Vg(^)(J|-). (25) 

A quadrupole buckling mode radiates gravitational waves at close to twice the angular 
frequency of torus (van Putten 2002a). These emissions are relatively powerful, repre- 
senting about 10% of the rotational energy of the black hole for the lifetime of the sys- 
tem. The associated mass inhomogeneity 5Mt in the torus assumes a value commensurate 
with the inferred luminosity in gravitational radiation. For quadrupole emissions, we have 
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Lgw ^ {?>2/b){MH / Rfi^Mr/ MhY , where uj ~ M^I^R-^/^ denotes the orbital angular fre- 
quency at a radius R, assuming approximately circular motion. The estimated gravitational 
wave-luminosity is hereby produced by SM^ — 0.5%Mh(R/5MhY^*. This corresponds to a 
mass inhomogeneity of 20% for a torus Mt — 0.2Mq around a black hole of mass Mh — 7Mq. 

Torus winds. The energy output in torus winds is a factor rj less than that in gravitational 
radiation, or 

B„.6xlO-erg(^)=(^). (26) 

These powerful torus winds may produce hypernova remnants in the host environment, e.g., 
a shell associated with a molecular cloud. They may be baryon loaded and deposit some 
torus matter onto the companion star as in the Brown et al. (2000) association of hypernovae 
to soft X-ray transients, and they may be important in coUimating baryon poor jets produced 
by the black hole. Finally, these winds arc potentially relevant in r-processes (Levinson & 
Eichler 1993). This suggests considering observational methods to determine E^^,, from which 
to determine the system parameter r/, and hence the frequency of gravitational radiation. 
Of potential interest are model dependent estimates of Ew from their role in collimating 
outflows, and calorimetry on hypernova remnants. The first can be pursued using existing 
data on GRB beaming angles, which suggests a value of 77 ~ 0.1. The second method is 
potentially more reliable, but awaits further study on the identification and observational 
aspects of hypernova remnants. 

Torus temperature. The energy output in thermal and neutrino emissions is a factor 5 less 
than that in gravitational radiation, or 

E, = 10-erg m ( ^] ( ^] . (27) 

Here, we refer to a fiducial value of 77 ~ 0.1 as before, as well as a value 5 ~ 0.15 or less. The 
latter is suggested by quadrupole radiation of gravitational waves which requires 5 < 0.16 at 
the threshold of Rayleigh stability, according to (22). This dissipation rate corresponds to 
a temperature of a few MeV (see eq. [29] below), which thereby produces baryonic winds. 
The torus winds considered here, therefore, are baryon rich. 



6. Pressure driven mass ejection from the torus 



It has been argued in §5 above that a fraction Ed/ Erot of the black hole-spin energy 
is dissipated in the torus, thereby heating it to a temperature in excess of a few MeV. 
This drives a powerful wind by the pressure gradients in the surface layers of the torus. 
The resulting mass ejection opens magnetic field lines that pass through the outer Alfven 
point, and folds some of those in the outer layers of the inner magnetosphere. This creates a 
change in poloidal topology in the form of a coaxial structure of open fiux tubes with opposite 
magnetic orientation (towards infinity). Because the torus is rotating and magnetized, the 
ejection of the wind is partially anisotropic. Specifically, mass flux is generally suppressed 
along magnetic field lines that are inclined toward the rotation axis, and enhanced along field 
lines that are strongly inclined away from the axis (Blandford & Payne 1982; Romanova et 
al. 1997), owing to centrifugal forces. The details of the outfiow depend on the heating and 
coohng rate of the corona, and on its structure, and analysis thereof is beyond the scope 
of the present discussion. In what follows, we provide a rough estimate for the mass flux 
expelled in the vertical direction, speculate on the implications for opening of magnetic fleld 
lines, and discuss the consequences of mass ejection in the equatorial plane for the energy 
extraction process. 



6.1. An estimate of the vertical mass flux 

A torus having a temperature in excess of a few MeV, and an average density ~ 
1Qi^'^{Mt/Q.IMq)/V2i gr cm~^, where V = lQ'^^V2i cm^ is the volume of the torus, cools 
predominantly by neutrino emission through electron and positron capture on nucleons. 
The cooling rate is given by (e.g., Bethe & Wilson 1985) 

eeap ^ lO'^iiTfo ergss-^cm-^ (28) 

where Tio is the average temperature in units of 10^° K. For a total energy dissipation rate 
of L = 10^^L52 ergs, we then find an average temperature of 

Tio ~ 2LJ/'(Mt/0.1Mo)-^/^ (29) 

As will be shown below, the mass fiux from the surface depends on the temperature 
profile in the neighborhood of the fiow critical point, where the density is well below the 
average. The latter can be determined in principle by equating the local heating and cooling 
rates, provided that adiabatic cooling there can be neglected (which we find to be justified 
only if the temperature exceeds ~ 2 x 10^° K). While the fraction of the black hole spin- 
energy dissipated in the surface layers is not well constrained, a fraction of the energy of 
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neutrinos escaping from the dense regions - well beneath the surface - will be deposited in 
the surface layers of the torus. The dominant absorption processes for the latter are neutrino 
capture on neutrons and protons and pair neutrino annihilation into electron-positron pair. 
These processes typically dominate at lower densities and higher temperatures. In spherical 
geometry, the heating rate is dominated by neutrino annihilation at densities below p = 
6 X 10^(L52/-R^g)^/^(-Ri//r)^ gr cm~^, where Ri, = 10^Ri,q cm is the radius of the neutrino 
production region (Levinson & Eichler 1993). The cooling rate due to electron-positron 
annihilation into neutrinos depends sensitively on temperature {e,^p oc T^), and dominates 
over the energy loss rate by electron and positron capture on nucleons at densities below 
p — 5 X lO^T^Q gr cm~^ (Levinson & Eichler 1993). Taking into account both processes, we 
find, up to a geometrical factor, that the critical point (see eq. [37] below) will be maintained 
roughly at the average temperature given by eq. (29). 

Now, the torus material should be a mixture of baryons and a light fluid (photons and 
electron-positron pairs in equilibriiim) . The light and baryonic fluids will be tightly coupled, 
owing to the large Thomson depth. Deep beneath its surface, the torus is in a hydrostatic 
equilibrium where the vertical gravitational force exerted on it by the black hole is supported 
by the baryon pressure Ph = rihkT. At baryon densities 

p,<^c.l0^^rfo grcm-^ (30) 

the light fluid pressure, 

= 2 X lO^^Ti^o dyncm-2, (31) 
exceeds the pressure contributed by the baryons. At sufficiently shallow layers, the hght 
fluid pressure gradient overcomes the vertical gravitational force, and the matter starts to 
accelerate. The transonic flow should pass through a critical point. To estimate the mass 
flux we calculate below the wind density and velocity at the critical point. 

To simplify the analysis we neglect general relativistic effects. Since we are merely 
interested in an estimate for the mass loss rate, we need only to consider the properties of 
the flow in the neighborhood of the critical point. Since the latter is located well within 
the light cylinder, we neglect rotation of the torus and the toroidal magnetic field. (This is 
not valid in general, but is used only in regards to mass ejection from the upper and lower 
faces of the torus.) Below, we find that the flow becomes mildly relativistic at the critical 
point. The MHD limit applies, in view of high electric conductivity in the coupled light plus 
baryonic fluids, whereby the streamlines of the flow are along magnetic flux-surfaces. Baryon 
number conservation, viz., V{nf,Up) — 0, where Up denotes the poloidal 4-velocity ans rib the 
baryon number density, and Maxwell's equation, V x E = 0, imply that the flux of baryons 
per unit magnetic flux is conserved: 

(nhUp/Bp)' = 0. (32) 



-26- 



Here Bp is the poloidal magnetic field, and ' denotes derivative along streamlines (i.e., V). 
In the limit of weak gravitational field, the projection of the momentum equation, T**^ = 0, 
on the poloidal direction, and the use of eq. (32) gives (Camenzind 1986; Takahashi et al. 
1990) 

Up{l - M-')u'^ = --t_[al{\nBpy + V''], (33) 

where 7^ = 1 + is the Lorentz factor of the flow, ip = GMn/rc^ is the gravitational 
potential, 

is the sound speed of the mixed fluid (measured in units of c), and M — Up/cg, with Cs — 
ttg/ ^yl — al being the sound four speed, is the corresponding Mach number. It is seen that 
the flow has a critical point at M = 1. We note that under the assumption made above, that 
the toroidal magnetic field can be ignored, this critical point coincides essentially with the 
fast magnetosonic point. We now suppose that the flow passes through this critical point, 
which should be true in the case of a wind expelled along open magnetic filed lines. There 
the right-hand side of eq. (33) must also vanish. The latter condition can be solved to yield 
the sound speed at the critical point: 



2 ^ ^' ^ GMh r' 
""'^ On Bp)' d'rl {luBpY ^''^> 

where Tc ~ a is the radius of the critical point. It is seen that the critical sound speed 
depends on the geometry of the field lines in the vicinity of the critical point. To obtain an 
order of magnitude estimate, let us assume that r' ~ Cc/tc, and (In Bp)' ~ where < Tc 
denotes the distance from the torus midplane to the critical point along streamlines. Then 



The critical density can be obtained now by employing eqs. (29), (34) and (36): 



The associated mass fiux is given by 

-1/2 / \ -2/3 



1 ^^^n ( \ f \ ' 1/2/^ , x-lr2/3 -2 -1 

[TOM-J [OAM^J ^-(^'^/^^^ S^^" ^ ' 



(38) 
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which corresponds to a mass loss rate of M ~ 1 x 10^° gr s~^ for a surface area oi A — 
O.la^ = lO^^Oy cm^. We conclude that during the ~ 30 s suspended accretion state the torus 
will be partially evaporated. 

Finally, we note that the Alfven Mach number at the critical point is 

Ma = (1^) Om{MT/0.1Mer'/'Lli'B;,\. (39) 

We thus conclude that, for our choice of torus parameters, the outflow is sub-Alfvenic (but 
not highly so) at the critical point. 

6.2. Creation of open field-lines 

The ejection of matter from the hot torus corona will result in the opening of some 
magnetic field lines in the outer layers of the inner torus magnetosphere (Fig. 10). By eq. 
(39), plasma which streams along field lines that extend to large distances quickly reach the 
Alfven point. These field lines become nearly radial several scale height above the torus, 
forming an open flux tube. For field lines that converge toward the rotation axis of the black 
hole the above analysis is probably inapplicable, as the centrifugal force cannot be ignored. 
Large pressure gradients in the corona would tend to push matter along some of the magnetic 
field lines in the outer layers of the inner torus magnetosphere. A combination of buoyancy 
and centrifugal forces may subsequently give rise to a twist of these field lines, some of which 
may ultimately fold and open, to form a region of oppositely directed magnetic field lines. 
Field lines thus created near the axis constitute the inner fiux tube that extends from the 
horizon to infinity: those anchored to the torus now extend to infinity and constitute an 
outer flux tube. The inner and outer flux tubes have opposite magnetic orientation and are 
separated by a charge and current sheet, whenever the outer tube carries a (super- Alfvenic) 
wind to inflnity. Reconnection may occur in the boundary layer, which is of interest in the 
rearrangement of the magnetosphere near the axis. This coaxial structure of open flux tubes 
is discussed in greater detail in the next section. 

A crucial issue is the opening angle of the inner flux tube on the horizon of the black 
hole, as it forms an artery for its spin-energy. We consider it probable that the opening 
angle depends only on the overall geometry of the system (black hole mass and radius of 
the torus) . A universal half-opening angle on the horizon gives rise to a standard fraction of 
black hole-spin energy released through the inner flux-tube. 
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Fig. 10. — The MeV temperature of the torus produces powerful torus winds. We propose 
that these winds open an outer layer of flux- surfaces, supported by the inner face of the 
torus. This is schematically indicated by a stretch, fold and cut of the scparatrix (dashed 
lines) between the inner and outer flux-surfaces. This leaves an open, inner tube of magnetic 
flux supported by the black hole and an open, outer flux tube of magnetic flux supported 
by the inner face of the torus. The inner tube is baryon-poor; the outer tube is baryon-rich. 
The lower /upper sections of this structure contain parallel/ antiparallel orientation of the 
poloidal magnetic field. The inner tube serves as an artery for baryon-poor outflow from the 
black hole. 
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6.3. Implications for energy extraction from the black hole 

The rate at which the torus catalyzes black hole spin-energy has been calculated in 
§5 under the assumption that the torus magnetosphere is force-free. However, ejection 
of appreciable mass flux along magnetic field lines that thread the horizon may alter the 
extraction process. Detailed study of ideal MHD flow in Kerr geometry is provided by 
Takahashi et al. (1990), and Hirotani et al. (1992). They analyze the conditions under 
which energy extraction occurs in a flow that starts with a zero poloidal velocity at the 
plasma source and is pulled inward by the black hole. They show that the MHD flow can 
carry a negative energy flux into the horizon (which is equivalent to the condition that the 
energy given by eq. [A9] is negative), provided that i) the angular velocity of magnetic field 
lines lies in the range < Qp < Qh, where Qh is the angular velocity of the black hole 
(which is identical to the condition found by Blandford & Znajek (1977)), and ii) the Alfven 
point is located inside the ergosphere, which implies that inflow of negative electromagnetic 
energy exceeds inflow of positive kinetic energy. The latter condition restricts the mass 
flux expelled along fleld lines connecting the torus and the horizon. The vertical mass 
flux estimated above assumes maximum extraction efficiency. We argue that inward mass 
ejection in the equatorial plane may be suppressed by the centrifugal barrier. The details 
are complicated by virtue of general relativistic effects. Whether this is sufficient to allow 
energy extraction is not clear at present. If not, it would mean that the rate at which the 
spin down energy of the hole is dissipated in the torus and, hence, its temperature, must be 
regulated by mass ejection. 



7. Structure of the inner flux-tube 

In the preceding section we argued that mass ejection from the torus opens magnetic field 
lines on the outer layers of the inner torus magnetosphere. This creates an open magnetic 
flux tube that extends from the horizon to infinity, surrounded by an outer fiux-tube that is 
anchored to the torus. The upper sections of this structure results from a fold, stretch and cut 
in poloidal topology, whereby the inner and outer tubes assume mutually antiparallel poloidal 
magnetic fleld. The lower section remains unfolded, leaving parallel poloidal magnetic flelds 
in the inner and outer flux-tubes (see Figs. 4, 10). This additional structure was not analyzed 
in (van Putten & Levinson 2002). In what follows, we describe the structure of the inner 
tube in some detail. 
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7.1. Asymptotic boundary conditions on the horizon and at infinity 

The inner flux tube satisfles slip/slip boundary conditions both on the horizon and at 
inflnity. It is well known that in the limit of inflnite conductivity, every magnetic surface 
must rotate rigidly (although the angular velocity of different flux surfaces should not be 
the same). Finite resistivity effects, however, may give rise to a differential rotation along 
magnetic flux tubes (implying E • B 7^ 0) . To make our analysis general, we shall allow for 
such a differential rotation, and denote the angular velocities of a given flux surface on the 
horizon and at inflnity by ftF+ and JIf-, respectively. 

As stated above, we anticipate the horizon half-opening angle 9h of the inner flux tube 
to be sufficiently small so that the small-angle approximation applies in solving Maxwell's 
equations. Prom equation (A20) we obtain 

27r J fV^de = A sin OFre = ^Bt- (40) 

It can be readily shown (e.g., Blandford & Znajek (1977)) that in the force-free hmit, viz., 
F/j.uj'^ — 0, the Boyer-Lindquist toroidal magnetic field, Bt, is conserved along magnetic fiux 
surfaces \E'(r, 9), and that the current contained within a magnetic fiux surface is given exactly 
by eq. (40), viz., /(^') = (1/2)5^ and, hence, is also conserved. Current conservation along 
streamlines is not guaranteed in general, however, even in the limit of infinite conductivity, 
since incrtial effects may give rise to cross field currents. At any rate, beyond the fast mag- 
netosonic point of the infiow (outfiow), the poloidal current becomes radial asymptotically, 
as it is carried purely by the infiowing (outfiowing) Goldreich- Julian charges (see below) . By 
employing the asymptotically frozen-in condition on the horizon, = {T,/ A){flH — ^F+)F^e 
(see eq. [A18] with /3 = —^h), we obtain the poloidal current fiowing through the horizon: 

Assuming the ZAMO radial magnetic field, B^ = —F^q/(T, sin 9), to be independent of 9 near 
the axis, we can calculate the total magnetic fiux through the horizon in one hemisphere: 
^ = 27rAtf, = J Br^yg^4,geed9d(j) ~ TrEsin^ 9Br — — Trsin^F^^. Eq. (41) then yields to leading 
order, 

Ih ^ -i^H - nF+)A^. (42) 

Likewise, the frozen-in condition at infinity gives Fre — ^.p-Ffj^e, leading in the small angle 
approximation to, 

7oo = (1/2) sin 9F^0nF- ^ -^f-A^. (43) 

The electric charge distribution in the inner fiux tube can be obtained from Maxwell's equa- 
tion: F*'^^ = 47rj*. At small angles this equation is given, to a good approximation, by eq. 
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(A19) in the Appendix (cf. Van Putten 2001). Assuming as before that the ZAMO radial 
magnetic field is independent of 9 near the axis, we obtain from eq. (A19) the Goldreich- 
Julian charge density: 



p^^^2.t^_iSlF + m^ose^ (44) 

where denotes the local angular velocity of the flux-surface at hand. Evidently, the electric 
charge changes along the inner flux tube from p^H — {^h — ^F+)Br/2'K on the horizon, to 
Peoo — —^F-Br/'^TT at inflnity, and vanish at the radius at which —(3 — Qp, corresponding 
to a locally zero angular momentum state of the flux-surface. This is illustrated in fig. 4. 
Using eq. (A20) and (44), one finds j'' = where f = 1 at infinity, and if = —cP' on 
the horizon (see Appendix). This implies that the current on the horizon and at infinity is 
purely due to convection of Goldreich- Julian charges. A similar conclusion was drawn by 
Punsly and Coroniti (1990). 



7.2. Two steady state limits 

Current closure of loo through the inner tube over the outer tube with equal magnetic 
flux of opposite sign gives rise to a differentially rotating inner tube, described by 

Qf+ - ^F- ^^H- 2VtT (45) 

(van Putten & Levinson 2002). This assumes the force-free limit, whereby current is con- 
served along flux-surfaces and Ih = loo- 

Alternatively, current closure of loo through the inner tube over the outer tube and 
the outer torus magnetosphere allows for the limit of a uniformly rotating inner tube. This 
corresponds to the limit of infinity conductivity. Combined with the force- free limit, this 
approach was used by Blandford & Znajek (1977) and Phinney (1983) to construct their 
force-free solutions and determine the efficiency of the energy extraction process. Near the 
axis Eqs. (42) and (43) yield to leading order, 

VLf+ = VLf- = VLh/2. (46) 

It will be appreciated that this implies nearly maximal energy extraction rate on the inner 
tube. 

Punsly & Coroniti (PC; 1990) argue that the assumption made by Blandford & Znajek 
(1977) and Phinney (1983), namely that the magnetosphere is force-free in the entire region 
between the horizon and infinity, violates the principle of MHD causality. Their argument 
relies on the observation that the inflow must become superfast on the horizon and, therefore. 
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cannot communicate with the plasma source region (e.g., the gap in the Blandford-Znajek 
model). They concluded that the use of the Znajek frozen- in condition on the horizon to 
determine Qp is unphysical, and that Qp must be determined by the dissipativc process that 
leads to ejection of plasma on magnetic field lines. This point is of interest, as the black hole 
has a finite, though secular lifetime of rapid spin as an inner engine to gamma-ray bursts. 
In what follows we reexamine this issue. 

We remark that MHD causality prohibits infinitely rigid structures (the Alfven velocity 
is bounded by the velocity of light). This is true in particular over distances much larger 
than the system size. Nevertheless, wc may examine a steady-state force-free limit Qf+ — 
QiT'- ~ ^h/'^ as a close approximation on scales comparable to the system size, modified 
only by a differentially rotating gap which injects the associated electric current in a region 
where frame-dragging /3 is pronounced. 

The current flow in the inner tube is created in a slightly differentially rotating gap 
in a neighborhood of = — /3 between two Alfven surfaces. The gap size (and hence 
current output) is determined by the local degree of differential frame-dragging in /3. As the 
magnetosphere around the black hole is transparent to /3, the gap remains in causal contact 
with the angular velocity of the black hole: a change in the asymptotic value —(3 — D,h on 
the horizon is associated with a change in f3 throughout the surrounding spacetime. This 
indicates that the gap is subject to changes Sfln under general conditions. 

If Ih 7^ -^oo) the inner tube develops a response by charge conservation: time- variable 
or in steady state, depending on the absence or presence of cross-field currents (i.e., je) 
between the inner and the outer tube. Conceivably, the plasma injection process maintains 
the required cross-field currents, and there exist multiple current closure paths. In view 
of (42) and (43), note that this implies additional differential rotation in the inner tube 
with accompanying dissipation processes. In the absence of cross-currents, there results a 
time-variable adjustment of the two Alfven surfaces which delimit the gap by application 
of Gauss' theorem to, respectively, the black hole plus lower section and the upper section. 
This causes adjustment of the current injected into the lower and the upper sections. In 
steady state, this recovers Ih — loo- This argument leaves open the possibility that the gap 
is time-variable, however, especially on the light crossing time-scale. 

7.3. Output power through the inner tube 

The inner tube releases black hole-spin energy at a certain rate, set by the horizon half- 
opening angle 6h shown in Figs. 4 and 10. Forementioned force-free hmits (45-46) define 



-33- 



bounds for the bipolar outflows of black hole-spin energy through the inner tube, i.e., 

Qj^+inH - nF+)Al <S^< \n],Al, (47) 

where 27r/lj„ is the net magnetic flux on the horizon associated with the inner tube. Hence, 
we have a fraction 

^ ^ J^H (48) 

of the rotational of the black hole (see eqs. [1] and [2]). 

The observed true GRB energies cluster around 3 x 10^° erg (Frail et al. 2001). This 
corresponds to Ej — 2 x lO^^erg x(0.16/e), where e denotes the efficiency of conversion 
of kinetic energy-to-gamma rays. In our model of GRBs from rotating black holes, this 
observational constraint introduces the small parameter 

, 10- (°-^] . (49) 



Erot \ M ^ ^ 

This corresponds to a horizon half-opening angle Oh > ^5°; 9h — 35° associated with the 
output from the gap in (45). 

We propose to identify 9h with the curvature in poloidal topology of the inner torus 
magnetosphere (van Putten 2002b), 6h ~ M/R within a factor close to one. In this event, 
we have 

Ei 1 fMV , , 



Erot 4 \ i? 

A spread in torus radius by a factor of about two corresponds to a spread in energies in 
baryon-poor outflows by about one order of magnitude, consistent with the observed spread 
in true GRB energies. 

The fate of the Poynting flux-outflow is determined by speciflc mechanisms for dissipa- 
tion, as these may arise downstream of the upper section of the open flux-tube. 



8. Structure of the outer tube 

The open held lines that are anchored to the torus form the outer flux tube. The 
baryon-rich material ejected from the torus, derived from the spin-energy of the black hole 
(see §6), flows along those open held hues to inflnity. We have estimated this to constitute 
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a substantial mass loss. In the ideal MHD approximation the angular velocity of the outer 
flux-tube is conserved along magnetic field lines (on the scale of the system) and, therefore, 
equals that of the torus, viz., Qp = Q-r- Above the Alfven point the ratio of toroidal and 
poloidal field in the outflow is (see eq [All] with = fix and v^j, r~^) 

= — ■ (51) 

As argued above, the direction of the poloidal magnetic held in the outer flux tube is opposite 
to its direction in the inner tube. Thus, if the torus and the black hole rotate in the same 
direction, as envisioned here, then the toroidal magnetic fleld F^e have opposite orientations 
in the inner and outer tubes. The Goldreich- Julian charge density in the outer tube is 

Pe = ' ^^^^ 

and is opposite in sign to the outflowing charges in the inner flux tube. These charges 
carry the current flowing along the open fleld lines from the torus. This region of the torus 
magnetosphere is equivalent to a pulsar magnetosphere. 

Wc conclude that the outer tube forms out of the outer layers of the inner torus magne- 
tosphere (the magnetic wall around the black hole), and joins the open fleld-lines in the outer 
torus magnetosphere in creating a super- Alfvenic baryon-rich wind to inflnity. It conceivably 
contributes to coUimation of the baryon-poor outflows in the inner flux-tube. 



9. The interface between the inner and the outer flux-tube 

The creation of the open magnetic flux tubes from the closed torus magnetosphere 
topologically represents folding of magnetic field lines in the upper section of the inner/outer 
fiux-tube, which gives rise to an antiparallel orientation of the poloidal magnetic field. This 
is accompanied by a cylindrical current and charge sheet that accounts for the jump in the 
electric and magnetic fields across the interface. The lower section of the inner/outer flux- 
tube which connects to the horizon has, instead, a parallel orientation between the poloidal 
magnetic fleld. In the perfect MHD limit, the properties of the interface are described by 
jump conditions as follow from Maxwell's equations, F^^.^ — 47rj*, i.e.: 

Anae = [nrsinOBr] = rsin^QT-B,.+ - rsine(Qjf/2)S^_, (53) 

where B^.^ (Br-) denotes the radial magnetic fleld near the interface in the outer (inner) 
flux tube. The poloidal and toroidal interface currents are, likewise. 
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The poloidal current (54) results beyond the Alfven point, where the wind transports angular 
momentum outwards to infinity. The onflow in the inner tube is expected to be relativistic 
{v'^ = !)• If the baryon rich outflow from the torus also becomes relativistic, then the latter 
equation imply that = cTg, namely the poloidal current in the boundary layer is solely due 
to the outflowing surface charges. The poloidal current sheet - possibly further carrying a 
poloidal current - is a potential site for reconnection of magnetic field lines, which would 
convert magnetic enegy in the inner flux tube into kinetic energy. 



9.1. Dissipation in the folded upper section of the flux-tubes 

The upper interface between the inner and the outer flux tubes is subject to antiparallel 
magnetic flclds and strong differential rotation (whenever Qt 7^ ^///2. It is therefore a 
potential site for reconnection of magnetic fleld lines, and may give rise to conversion of a 
fraction of the magnetic enegy in the inner flux tube into kinetic energy. In order to asses the 
fraction of magnetic energy which is converted in the interface, a detailed reconnection model 
is required. Here, we point out that the reconnection time is limited by the crossing time 
af an Alfven wave across the inner tube, which is typically very short. If this limit applies 
then an appreciable fraction of magnetic energy will be dissipated. (The reconnection rate 
could be much smaller.) In any case, the power extracted from the hole by the inner tube 
is uncertain, because the angular velocity of the inner tube on the horizon is unknown. In 
or near the lowest energy state of the gap, the inner tube mediates a power S^n, as given by 
eq. (47). 

The electric flelds produced by the reconnection process in the interface may inject 
plasma, including baryons that should be present in the boundary layer, into the inner tube. 
This would lead to mass loading of the inner tube and additional conversion of the Poynting 
flux. The amount of baryonic contamination by this process is yet an open issue. 

Additional dissipation may arise from a differential rotation of the inner flux tube, as 
discussed in (van Putten & Levinson 2002). It is clear that if the system fluctuates over a 
time scale At, then nonlinear disturbances may induce differential rotation over length scales 
< cAt (which corresponds to 10^ gravitational radii for the entire life time of the system). 
Such electromagnetic fronts should be highly dissipative by virtue of the large parallel electric 
flelds associated with the differential rotation of the magnetic flux tubes. It is anticipated 
that, under the conditions envisioned here, copious pair and photon productions would ensue 
inside the differentially rotating fronts. 
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10. Observational consequences 

We have described in some detail the topology, lifetime and emissions of black hole- 
torus sytems in the suspended accretion state. Fig. 11 summarizes the energy transport and 
conversion of black hole-spin energy, catalyzed by the torus. Most energy is dissipated in 
the horizon, while the major output is converted by the torus into gravitational radiation, 
winds and neutrino emissions (as well as thermal emissions). A minor output is released in 
baryon-poor outflows as input to the observed GRB-afterglows. The analysis is based largely 
on equivalence to pulsar magnetospheres. The suspended accretion state develops against a 
magnetic wall around rapidly rotating black holes. This is based on a uniform magnetization 
of the torus, represented by two oppositely oriented current rings in a torus formed in black 
hole-neutron star coalescence and core-collapse hypernovae. 

The proposed emissions in gravitational radiation are powered by the spin-energy of 
the black hole, which renders these emissions a powerful new source for the upcoming grav- 
itational wave-experiments. This output is contemporaneous with GRBs from baryon poor 
outflows, and hence have the same durations as the redshift corrected durations of long 
GRBs - tens of seconds, due to the condition of magnetic stability of the torus. 

Calorimetry on the emissions in gravitational radiation provides a rigorous compactness 
test for Kerr black holes (van Putten 2001a), which can be pursued by upcoming gravita- 
tional wave-experiments. We emphasize that our predictions on the energies and duration 
of emissions in gravitational radiation are robust and are independent of an association with 
GRBs. Hypernovae could be over-abundant as burst sources of gravitational radiation, with 
only a small fraction making successful GRBs. Calorimetry on wind energies provides a 
method for constraining the angular velocity of the torus and, hence, its quadrupole emis- 
sions in gravitational radiation. Possible channels are recently observed X-ray line emissions 
in GRB 011211 (Reeves et al. 2002), hypernova remnants - which remain to be identified 
- and by association with soft X-ray transients with chemically enhanced companion stars 
(Brown et al. 2000). 

X-ray and radio shells that may result from the interaction of the baryon rich torus 
winds with ambient matter, as well as the gamma-ray emission expected to be produced by 
the baryon poor outflows in the inner tube, should be sought for (Levinson et al. 2002). 

Below we summarize our findings and draw additional conclusions: 

1. Most of the black hole- luminosity is incident on the torus by topological equivalence 
to pulsar magnetospheres, when the black hole spins rapidly. We identify a new magnetic 
instability, which defines a minimum lifetime T of tens of seconds of the spin of the black 
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Fig. 11. — Overview of energy transport and conversion of black hole-spin energy catalyzed 
by a uniformly magnetized torus. Most of the spin-energy is dissipated in the horizon - an 
unobservable sink of energy - at a rate which is limited by a new magnetic stability criterion 
for the torus. This process effectively defines the lifetime of rapid spin of the black hole. 
Most of the spin-energy released is incident on the inner face of the torus, while a minor 
fraction forms baryon-poor outflows through the inner flux-tube to infinity. We associate 
the latter with the input to the observed GRB-afterglow emissions. The torus converts its 
input primarily into gravitational radiation and, to a lesser degree, into winds, thermal and 
neutrino emissions. Direct measurement of the energy and frequency emitted in gravitational 
radiation by the upcoming gravitational wave-experiments provides a calorimetric compact- 
ness test for Kerr black holes (dark connections). Channels for calorimetry on the torus 
winds are indicated below the dashed line, which are incomplete or unknown. They provide 
in principle a method for constraining the angular velocity of the torus and its frequency of 
gravitational radiation. This is exemplified by tracing back between torus winds and their 
remnants (dark connections). As the remnant envelope expands, it reaches optical depth of 
unity and releases the accumulated radiation from within. This continuum emission may 
account for the excitation of X-ray line emissions seen in GRB 011211, which indicates a 
torus wind energy of a few times lO^^erg. Matter ejecta ultimately leave remnants in the host 
molecular cloud, which remain to be identified. Torus winds may further deposit a fraction 
of their mass onto the companion star (Brown et al. 2000), thereby providing a chemical 
enrichment in a remnant soft X-ray transient. 
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hole, corresponding to a large system parameter T/M ~ 10^. 

2. At MeV temperatures of the torus, a torus with a minor radius of less than about 
the linear size of the black hole becomes unstable to the formation of a finite multipole mass 
moments by the Papaloizou-Pringle instability and, conceivably, other modes. The torus 
hereby converts a major fraction Egw/Erot — 10% into gravitational radiation, conceivably 
at multiple lines. At a source distance of lOOMpc, these emissions over N = 2 x 10^ periods 
give rise to a characteristic strain amplitude VNhchar — 6 x 10~^^ (Coward et al. 2002). 

3. A minor fraction Ej/Ej-ot ~ 10~^ of black hole-spin energy is released in baryon poor 
outflows. We attribute this fraction to curvature in poloidal toplogy, whereby Ej/Erot — 
(1/4) (M/i?)'^ is standard for a universal horizon half-opening angle 6h — M/R of the asso- 
ciated open magnetic fiux-tube. These outflows are probably not high-sigma - the ratio of 
Poynting flux-to-kinetic energy flux - owing, in part, to magnetic reconnection in an interface 
with baryon-rich winds flowing along the surrounding outer flux-tube. 

4. Causality in the process of spin-up of the torus by the black hole follows by topological 
equivalence to pulsar magnetospheres (van Putten 1999). In response, the black hole spins 
down by conservation of energy and angular momentum. This establishes causality in the 
extraction of black hole-spin energy by horizon Maxwell stresses as proposed by Blandford 
& Znajek (1977). Causality in the formation of baryon poor outflows from the open mag- 
netic flux-tube is due to transparency of the magnetosphere to the gravitational field: the 
associated current injection is subject to current continuity between asymptotic boundary 
conditions on the horizon and at infinity by Gauss' theorem, and regulated by differential 
frame-dragging and current closure at infinity. 

5. Calorimetry on the predicted energy output in gravitational radiation by the upcom- 
ing gravitational wave experiments provides a method for identifying Kerr black holes as 
objects in nature when 27r /q^""" fgwdE > 0.005 (van Putten 2001a), where fgw denotes the 
quadrupole frequency of gravitational waves. Current experiments consist of layer interfero- 
metric detectors LIGO, VIRGO, TAMA and GEO, bar and sphere detectors. An individual 
source is band limited in gravitational radiation in terms of a frequency sweep of about 10%, 
corresponding to emission of the first 50% of the output in gravitational radiation from a 
maximally spinning Kerr black hole. Collectively, black hole-torus systems are conceivably 
luminous in multiple frequencies with broad distributions owing to a spread in black hole 
mass. 

6. Frequencies of quadrupole gravitational radiation fgyj ~ 470Hz {E^/A x lO^^erg)^/^ 
{7Mq/Mh)^^'^ can be predicted by calorimetry on the torus wind energies Eyj. Calorimetry 
on X-ray line emissions point towards frequencies around 500Hz (van Putten 2002b). 



-39- 



7. Calorimetry on Eyj may also be pursued by calorimetry on hypernova remnants, 
given the observed supernova association (e.g., Bloom et al. (2002)). Our model suggests 
the ejection of the remnant stellar envelope as a shell with kinetic energy 0.5f3Eyj ~ 3 x 
10^^ cTg{l3 / 0. 1) {Mh/ 10 Mr..,) [r]/ 0.1). Here, (3 denotes the initial radial velocity (3 relative to 
the velocity of light, in response to the impact of E^^ from within. For a supernova association 
to molecular clouds, see, e.g., Chu & MacLow (1990); Wang & Helfand (1991); a hypernova 
or gamma-ray burst association to molecular clouds has been considered by Efremov et 
al. (1998); Wang (1999); Dunne et al. (2001); Lai et al. (2001); Chen et al. (2002); Price 
et al. (2002). Determining kinetic energies of these remnants can be pursued analogously 
to studying supernova remnants, in the radio and X-ray. Chandra-X observations may 
hereby identify hypernova remnants in X-ray bright (super)shells around a black hole-binary, 
possibly in the form of a soft X-ray transient. 
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A. Ideal MHD in Kerr geometry 

We express the Kerr metric in Boyer-Lindquest coordinates with the following notation: 

ds^ = -a^dt^ + Co'^{d<j) + (3dtf + ^cir^ + p^dO'^, (Al) 

where a — pVA/T, is the lapse function, o;^ = (E^/p^) sin^ 9, —f3 — 2aMr/E^ is the angular 
velocity of a ZAMO with respect to a distant observer , with A = + — 2Mr, = 
+ cos^ 9, and = (r^ + a^)^ — a^A sin^ 9. The parameters M and a are the mass and 
angular momentum per unit mass of the black hole. 

We denote by n, p, p, h — {p+p)/n, the proper particle density, pressure, energy density, 
and specific enthalpy, respectively The stress-energy tensor then takes the form: 

T"^ = hnu'^u^ - pg''^ + -^(F""F^ + \g''^F^), (A2) 

Air 4 

where u"' is the four-velocity, and F^^ is the electromagnetic tensor which satisfies Maxwell 
equations. The dynamics of the MHD fiow is then governed by the equation 

T'^;^ = 0, (A3) 

subject to conservation of particle fiux, 

(nu^),^ = -^{^nu^),^ = 0. (A4) 

The ideal MHD assumption (i.e., infinite conductivity) imposes the additional constraint, 

u^'F^p = 0. (A5) 

The stationary axisymmetric flow considered here is characterized by two Killing vectors: 
= dt and x'^ = d^. By contracting with the stress-energy tensor we can construct the 
energy and angular momentum currents: = T^^^ = T", and = T^X^ = which 
are conserved, viz., 

/:% = = 0. (A6) 

Equations (A4), (A5), and (A6) together with the homogeneous Maxwell equations 
admit 4 quantities that are conserved along magnetic flux surfaces, the shape of which is 
given by a stream function ^(r, 9): The particle fiux per unit magnetic fiux, 

nm = ^ = (A7) 
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The angular velocity of magnetic field lines, 

= = Fte/Fe^. (A8) 

The total energy and angular momentum per particle carried by the MHD flow, 

Em = hut - ^^^fF'^, (A9) 
47r?7 

Li^) = -hu, - ^F^'. (AlO) 

The above equations also yield the relation 

^ = (All) 

where v"^ = u'f'/u^, and v'^ = /u^ are the corresponding components of the 3- velocity. 
Equations (A9) and (AlO) can be used to express -u*, v"^ , and F^^ in terms of L, and the 
Alfven Mach number, defined by = A'Khrf /n. One finds 

= ^ k^TJp ' ^^^^^ 

hu' = h{-u, + Pu,)/a' ^{E- QfL) - ^^l^^l^^^y (A13) 

4> ^ -c^QH^FjE - npL) + M\gt^E + guL) 

-a^g^4,{E-nFL)-M\g^^E + gt^L) ' ^ ' 

where = g^ + 2gt^VLp + g^^VPp. Finally, using the normalization condition, u^^u^j^ = — 1, 
we obtain for the poloidal velocity, defined as — UrU^ + uqu^ , the relation 

ul + l^ {au'f - {u^/ujf. (A15) 

Note that Up is the poloidal 4- velocity as measured by a ZAMO, and au^ is the corresponding 
t component of the 4- velocity. As seen, the Lorentz factor of the inflow in the ZAMO frame 
approaches V = au^ ~ a^^ on the event horizon. 

Consider first the behavior of the solution near the horizon. There a ^ and so 

, {^t^E + guL) 



{g^^E + gt^L) 



Hh- (A16) 



The poloidal velocity is radial on the horizon, implying Up — > grrU^u^ ■ Using eq. (A15) we 
find, 

^vT/u* ^ a/yfg;;^ -l^/{r'^ + a^). (A17) 
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Substituting the above results into eq. (All), we finally obtain, 

in + (3). 



Fr9 

Frhft 



(A18) 



Eq. (A18) gives the boundary condition on the horizon. 



The Boyer-Lindquist electric charge density in the flow is determined from Maxwell's 
equation: F*'*^ = 47rj*. It can be readily shown that if the poloidal component of the 
magnetic field is radial on the event horizon (which is the case in the force-free limit of an 
uncharged blck hole), then near the horizon the later equation reduces to. 



d_ 



sin 9 



a" 



= 47rA 
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(A19) 



The poloidal current follows from the equation, 

d / sin ^ _ 



de 



(A20) 



Using eqs. (All), (A16), (A19), and (A20), we finally recovers the result (Punsly & Coroniti 
1990)) 

f = (A21) 



where Pe — ^-nd Vp — Up/ui being the three velocity as measured by a ZAMO. 



